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Introduction 


Currently  breast  cancer  screening  is  performed  using  x-ray  mammography,  including 
both  the  screen-film  based  mammography  and  digital  mammography.  These  two  x-ray 
mammography  modalities  all  rely  on  small  differences  in  the  linear  attenuation  coefficients 
between  adipose  and  glandular  tissues  to  produce  the  mammogram.  Breast  is  one  of  the  most 
difficult  human  organs  for  x-ray  imaging,  since  it  consists  of  tissues  with  similar  densities  and 
attenuation  coefficients,  resulting  in  low  contrast.  However,  as  x-ray  wave  passes  through 
tissues,  x-ray  undergoes  phase-shifts  as  well.  The  amount  of  the  phase  shift  is  determined  by 
tissue’s  refractive  index  and  the  ray  path  length.  The  refractive  index/?  for  x-ray  is  complex  : 
n  =  \-S  +  ijB  ,  where 5  is  the  refractive  index  decrement  determined  by  x-ray  wavelength  and 
tissue  electron  density,  is  responsible  for  x-ray  phase  shift,  and  j3  is  the  imaginary  part  of  the 
refractive  index  and  it  is  responsible  for  x-ray  attenuation.  The  amount  of  x-ray  phase  shift  from 
a  piece  of  tissue  is  given  by  <f>  =  -{in  /X)  J  S(s)ds  =  -/Lre  |  pe(s)ds  ,  where  A  denotes  x-ray 

wavelength,  re  the  classical  radius  of  electron,  pe  tissue's  electron  density,  and  the  integral  is 
taken  over  the  ray  path.  Consider  the  infiltrating  ductal  carcinoma  (I DC)  in  breast,  we  found  that 
its  5  is  about  1 000  times  greater  than  its /?  for  x-rays  of  1 0  keV-1 50  keV  [30].  Similar 
conclusions  can  be  drawn  for  other  biological  tissues.  It  is  ironical  that  the  current 
mammography  techniques  are  designed  to  image  a  tissue's  p  but  not  its£  .  Hence  the  phase 
sensitive  x-ray  imaging  has  potential  for  greatly  increasing  x-ray  imaging  sensitivity  and 
specificity  and  reducing  radiation  doses  associated  with  the  imaging. 

The  approaches  for  phase-sensitive  x-ray  imaging  currently  under  investigation  fall  into 
three  broad  categories:  interferometric  imaging,  diffraction-enhanced  imaging,  and  inline  phase- 
sensitive  imaging.  The  technique  being  developed  in  this  project  falls  into  the  category  of  inline 
phase-sensitive  imaging,  which  is  in  general  based  on  the  inline  holography  principle  in  optics. 

In  inline  imaging  the  diffracted  x-ray  with  shifted  phases  forms  bright  and  dark  fringes  at  tissue 
boundaries.  The  edge  enhancement  thus  generated  relies  on  the  spatial  coherence  of  the  x-ray 
source,  the  optimal  positioning  of  detector,  and  the  x-ray  energy.  The  inline  phase-sensitive 
imaging  requires  neither  expensive  crystal  analyzers  nor  the  fabrication-challenging  gratings  as 
other  techniques  do.  For  the  inline  phase-sensitive  imaging,  the  spatial  coherence  of  x-rays  is 
more  important  than  its  temporal  coherence.  Furthermore,  the  spatial  coherence  offered  by  a 
synchrotron  radiation  source  can  also  be  obtained  with  a  x-ray  tube  with  a  small  focal  spot  size 
and  relatively  large  source-to-detector  distance.  Of  course  this  approach  for  achieving  x-ray 
spatial  coherence  is  constrained  by  the  requirement  of  sufficient  photon-flux  required  for 
achieving  short  exposure  time  for  clinical  imaging. 

The  long  term  objective  of  the  project  is  to  develop  a  low-dose  and  quantitative  phase  x- 
ray  imaging  technique  for  facilitating  breast  cancer  detection.  The  Specific  Aims  of  the  project 
of  the  three-year  period  are:  (1)  Develop  prototype  phase-imaging  system  enabling  the  phase 
retrieval,  that  is,  the  reconstruction  of  objects  phase-maps.  The  system  hardware  comprises  a 
micro-focus  tube  operating  at  high  tube  voltages,  a  high  resolution  photostimulable  phosphor 
plate  (CR-plate)  based  detector  system.  The  core  algorithms  for  breast  phase-map 
reconstruction  will  be  developed  to  retrieve  a  breast  phase  map  from  a  single  recorded  image. 
(2)  Validate  the  accuracy  of  the  reconstructed  tissue  projected  electron  densities;  validate  the 
many-fold  radiation  dose  reduction  achieved  with  the  proposed  system;  conduct  subjective 
measurements  to  characterize  the  performance  of  the  proposed  system. 
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In  the  first  year  of  this  project,  as  planned  in  the  Statement  of  Work,  we  performed  two 
tasks:  (A).  Develop  the  phase  retrieval  algorithms  for  future  phase  imaging  with  breast 
phantoms;  (B).  Design  and  build  the  system  hardware  for  future  phase  imaging  with  breast 
phantoms. 


(A).  Develop  the  phase  retrieval  algorithms 

Phase  contrast  x-ray  imaging  requires  a  x-ray  source  proving  highly  spatially  coherent 
illumination  of  the  imaged  object.  With  a  micro-focus  x-ray  tube  the  achievable  x-ray  spatial 
coherence  is  very  limited,  hence  the  phase-contrast  imaging  with  a  x-ray  tube  source  can 
generate  only  a  limited  phase-contrast  mixed  with  the  attenuation  contrast.  The  contrast 
enhancement  at  tissue  interfaces  and  boundaries  generated  by  the  dark-bright  diffraction 
fringes  are  still  limited.  On  the  other  hand,  the  phase  imaging  technique  that  we  proposed  goes 
behind  the  phase-contrast  imaging.  By  means  of  phase  retrieval  from  acquired  phase-contrast 
images,  one  can  disentangle  the  high  phase  contrast  from  the  mixed  contrast  mechanisms 
presented  in  the  phase-contrast  image.  Moreover,  a  retrieved  tissue  phase  map  provides  a 
map  of  imaged  object’s  projected  electron  densities  (the  ray  integrals  of  electron  densities  )  as 
well. 


The  ch  allenge  for  dev  eloping  br  east  x  -ray  phase  i  maging  is  how  t  o  overcome  the 
difficulty  f  or  i  mplementing  phase  r  etrieval  with  breast  i  maging.  In  or  dert  o  co  mpletely 
disentangle  tissue  phase-shifts  from  the  mixed  contrast  in  acquired  phase-contrast  images,  in 
general  one  needs  to  acquire  at  least  two  and  mostly  four  to  five  images  with  different  object- 
detector  distances.  This  requirement  is  imposed  by  the  x-ray  wave  nature  and  its  wave 
equation.  This  cumbersome  procedure  delivers  higher  radiation  dose  to  breast  and  makes  the 
images  su  sceptible  t  o  t  he  m  otion  a  rtifacts.  Therefore,  an  i  important  c  hallenge  i  s  to  find  a  n 
effective  and  low-dose  approach  for  phase  retrievals. 

Phase  retrieval  from  a  single  phase  contrast  image  is  the  key  strategy  for  this  project. 
We  previously  studied  how  high  energy  x-ray  wave  propagates  in  tissues.  Analyzing  the  tissue 
attenuation  data,  we  realized  that  the  soft  tissue  attenuation  cross  sections  are  approximate  by 
that  fr  om  x-ray  i  ncoherent  s  cattering  f  or  x-rays  of  abou  t  60  keV  or  higher.  Under  t  his 
circumstance,  we  found  that  soft-tissue  attenuation  and  soft-tissue  phase  are  all  related  to  the 
projected  electron  density,  we  call  this  new  notion  as  the  phase-attenuation  duality.  Based  on 
this  notion  of  the  phase-attenuation  duality,  we  discovered  a  new  innovative  way  to  retrieve  a 
phase-map  from  j  ust  a  si  ngle  recorded  image  [Wu  et  al.  2005].  In  this  way  t  he  pr  oblems 
associated  with  the  requirement  of  multiple-image  acquisition  for  phase-retrieval  is  eliminated. 
Recently  we  noted  that  our  previous  phase  retrieval  method,  though  performed  phase-retrieval 
very  well  in  general  from  only  a  si  ngle  phase -contrast  image,  is  not  accurate  enough  for  the 
cases  with  h  igh  sp  atial  resolution  and  large  o  bject-detector  di  stances.  I  n  order  to  develop 
robust  algorithms  with  high  accuracies  in  breast  phase  map  retrieval,  in  this  period  of  the  project 
we  investigated  and  developed  a  new  and  more  accurate  phase  retrieval  method  based  on  the 
phase-attenuation  duality. 

Method:  The  physics  behind  our  method  developed  in  this  project  can  be  explained  as 
follows.  The  phase  contrast  is  formed  by  the  phase-shifted  x-ray's  Fresnel  diffraction  from  the 
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object  to  the  detector.  The  enhanced  tissue  phase  contrast  is  manifested  as  the  formation  of 
bright  and  dark  diffraction  fringes  in  the  measured  images.  Our  task  is  to  determine  the  map  of 
tissue  phase-shifts  from  these  enhanced  image  intensity  variation.  We  noticed  that  in  the 
transverse  spatial-frequency  space  (the  transverse  k-space)  the  Fresnel  diffraction  can  be 
thought  as  a  x-ray  wave-front  filtering  with  the  Fresnel  propagator  exp(inhzu2) ,  where!  is  x- 
ray  wavelength  and  u  is  the  wave-front's  transverse  spatial  frequency  vector,  and  z  is  the 
object-detector-distance.  So  the  phase  of  the  propagator,  i.e.  n. Zzfi2 ,  determines  how  strongly 
the  exiting  wave-front  will  be  diffracted  in  the  propagation.  Our  previous  phase  retrieval  formula 
was  derived  by  linearizing  Fresnel  propagator  phases  for  simplifying  the  x-ray  wave  equation,  or 
the  transport-of-intensity  equation  [Wu  et  al.  2005],  Therefore,  in  order  to  develop  robust 
algorithms  with  high  accuracies  for  breast  phase  map  retrieval,  one  cannot  linearize  Fresnel 
propagators,  rather  one  should  consider  the  full  exponential  propagators  instead,  with  a  price  of 
dealing  a  much  more  complicated  x-ray  wave  equation.  We  then  used  the  phase-attenuation 
duality  to  delineate  the  complicated  equation  for  deriving  a  new  phase  retrieval  formula. 

Result:  Our  research  efforts  through  the  first  period  of  this  project  eventually  paid  off:  we 
successfully  developed  a  new  phase  retrieval  method  with  a  full  account  of  the  nonlinear 
Fresnel  propagators.  This  new  phase  retrieval  method  provides  the  scientific  foundation  and 
framework  for  developing  new  and  more  accurate  phase  retrieval  algorithms  for  this  project. 

The  computer  simulation  demonstrated  superiority  of  this  new  method  over  our  previous  duality- 
based  phase-retrieval  formula  in  reconstructing  the  high  frequency  components  of  imaged 
objects.  In  addition,  this  new  phase  retrieval  method  lays  the  foundation  for  implementing  high- 
resolution  phase-sensitive  volumetric  imaging  (CT  and  tomosynthesis)  of  soft  tissue  objects  as 
well.  Our  work  was  published  in  Optics  Express  in  June  2009  [Wu  et  al.  2009],  The  paper  is 
enclosed  in  the  Appendix  of  this  report  for  review. 

Emanating  from  this  result,  as  a  by-product  of  this  project,  we  have  applied  our  phase 
retrieval  method  to  the  phase  sensitive  tomography  as  well.  Previously  several  group  reported 
that  in  the  phase-contrast  CT  experiments  images  reconstructed  with  the  standard  filtered 
backprojection  algorithms  suffered  from  erroneous  artifacts  such  as  the  streaks  and  negative 
densities.  We  believed  that  these  artifacts  arise  from  the  fact  that  x-ray  does  not  propagate 
along  straight  lines  in  phase-contrast  projections,  because  x-ray  refracts  and  diffracts  at  tissue 
interface  and  boundaries.  In  order  to  correctly  deal  with  the  issue  of  phase-contrast 
tomography  reconstruction,  one  should  first  retrieve  the  phase  maps  for  all  projections,  then 
reconstruct  the  image-slices  of  tissue  attenuation  coefficients  and  electron  densities, 
respectively.  However,  multiple  image-acquisitions  required  for  phase  retrievals  make  the 
implementation  very  cumbersome,  since  multiple  phase-contrast  images  have  to  be  acquired  for 
each  projection  view  during  a  source  sweep.  Moreover,  radiation  dose  to  breast  will  be 
multiplied  for  the  multiple  image  acquisition  per  projection  view.  Therefore,  the  phase  retrieval 
method  reported  here  for  this  project  can  greatly  facilitate  the  artifact-free  phase-sensitive  CT. 

By  computer  simulation  we  showed  how  the  artifacts  get  generated  in  phase-contrast  CT  and 
how  our  phase  retrieval  method  developed  for  this  project  provided  artifacts-free  accurate  CT 
image  reconstructions.  We  reported  this  progress  at  the  2009  AAPM  annual  meeting  in  July 
2009.  The  peer-reviewed  abstract  of  this  oral  presentation  was  published  in  Medical  Physics 
and  is  enclosed  in  the  Appendix  as  well. 


(B).  Design  and  build  the  system  prototype 
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The  system  design  study  and  system  hardware  integration  have  been  conducted  in  this 

period. 

(1).  Study  on  the  system  design 

Determine  Imaging  Geometry  :  Different  from  conventional  x-ray  imaging,  the  phase 
contrast  imaging  needs  adequate  spatial  coherence  of  x-rays  illuminating  the  breast,  hence  it 


X-Ray  Source 


Detector 

3 


R1  R2 

Figure  1.  Geometry  Setting  of  the  phase  imaging  prototype.  The  source- 
object  distance  is  denoted  by  and  the  object -detector  distance  by  R2. 

requires  a  micro-focus  tube  as  the  source  ,  and  a  large  source-object  distance  Ri  as  is  shown  in 
the  Fig.  1 .  The  object  to  detector  distance  R2  should  be  sufficiently  large  as  well  to  allow 
phase-shifted  x-rays  to  diffract  as  forming  a  phase-contrast  image.  Considering  future  clinical 
application,  we  use  a  source  to  image-receptor  distance  (SID)  of  1 .5  m  -  1 .8  m.  Selecting 
optimal  values  of  the  distances  R t  and  R2  is  crucial  for  successful  phase  contrast  imaging.  This 
is  because  the  imaging-geometry  should  optimize  the  phase  contrast  visibility  by  balancing  the 
conflict  requirements  of  x-ray  spatial  coherence,  large  diffraction  fringes  and  finite  sizes  of  the 
focal  spot  size  and  finite  detector  pitch.  In  a  previous  work  we  found  that  the  phase-contrast 
that  could  be  imaged  is  proportional  to  the  relative  phase-visibility  factor  (RPF)  and  its  modulus 

is  \RPF(u)\  =  c2h2MTFGU  (u/M)MTFdet(u/M)^^R2u2 /ME2)s(E)dE^  ,  [Wu  etal.  2004], 

In  this  formula  h  is  the  Planck  constant  and  c  the  speed  of  light,  u  the  spatial  frequency  vector  of 
a  fine  object  structure  to  be  imaged,  M  is  the  magnification  factor,  E  x-ray  photon  energy,  and 
S(E)  is  the  x-ray  spectrum.  In  addition,  MTFGU  (u) denotes  the  MTF  of  the  geometric 

unsharpness  due  to  the  focal  spot  size,  and  MTFdet  (7)  is  the  MTF  of  the  detector.  For  a  given 
object  structure  to  be  imaged,  the  larger  |y?JP7^(77)|  a  system-configuration  achieves,  the  more  the 
phase-contrast  manifests  with  this  configuration,  though  the  actual  phase-contrast  exhibited  in 
an  image  depends  on  the  imaged  object’s  phase  variation  as  well.  The  higher  \RPF(u)\  is  the 
more  the  phase-contras  manifests.  Therefore  the  relative  phase-visibility  factor 
\RPF(ii)\  provides  a  very  useful  tool  for  designing  a  phase-contrast  imaging  system.  We 

studied  the  \RPF(u)\  for  this  project  and  a  typical  data  from  the  study  is  plotted  as  shown  in  Fig. 

2  below,  where  we  assumed  a  source-detector  distance  SID  =  1 .5  m,  a  uniform  circular  focal 
spot  of  50  in  diameter,  and  a  detector  pitch  of  43.75  |xm .  These  parameters  are  according  to 
the  specs  of  the  available  micro-focus  tube  and  imaging  detector.  In  addition,  we  assumed 
60.5  keV  x-ray  and  a  fine  object  structure  of  15  line-pairs/mm.  We  found  that  the  MTF  of  the 

geometric  unsharpness  \sMTFGU  (u)  =  \2Jl(7rf(M -l)|7|/M)]/[;r/(M -1)|m|/m] ,  where /  is 
the  diameter  of  the  focus  spot,  and  Jx  (x)  is  a  Bessel  function  of  the  first  kind.  We  used  the 
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sinc-function  for  modeling  the  MTFdet  («) .  Fig.  2  shows  the  computed  \RPF(u)\  plot  for  a  fine 

structure  of  15  line-pairs/mm.  We  found  that  the  \RPF(u)\  peaks  for  R-i  =  0.65  m,  and  R2=  0.85 

m,  that  is,  with  the  magnification  factor  M  =  2.31  for  60.5  keV  x-ray.  This  maximal  phase 
visibility  is  a  83%  increase  compared  to  the  phase  visibility  with  R!  =  1.1  m  and  R2=  0.4  m  (M  = 
1 .36.)  The  results  show  as  well  that  the  phase  visibility  will  decrease  if  the  magnification  factor 
exceeds  2.35.  Therefore,  our  reference  design  adopts  system  geometry  parameters  ranging 
from  Ri  =  0.7  m,  and  R2=  0.8  m  to  Ri  =  0.65  m,  and  R2=  0.85  m. 


RPF(u^)  vs.  Magnification 


Figure  2.  Relative  phase  contrast  visibility  \RPF(u)\  varies  for  different  imaging 

geometry  parameters  (Ri  and  R2 )  of  the  prototype  system.  In  this  plot  we  assumed  a 
source-detector  distance  of  1 .5  m,  a  uniform  circular  focal  spot  of  50  jam  in  diameter, 
and  a  detector  pitch  of  43.75  pm ,  60.5  keV  x-ray  and  a  fine  object  structure  of  15 
line-pairs/mm. 


According  to  above  analysis  and  the  hardware  availability  in  the  lab,  the  system 
hardware  components  to  be  employed  for  this  project  are  described  as  follows. 

X-Ray  Tube:  We  employ  a  micro-focus  tube  (Model  L8121-01,  Hamamatsu  Photonics) 
with  integrated  high  voltage  power  supply.  The  tube  has  a  tungsten-target  and  a  uniform 
circular  focal  spot  of  50  pm  in  diameter.  The  tube  can  operate  at  1 50  kVp  with  maximal  tube 
current  0.5  mA.  The  target  angle  is  25°.  Our  test  showed  that  with  such  a  tube  operating  at  120- 
150  kVp  breast  phantom  images  can  be  acquired  with  a  exposure  of  few  seconds. 

Image  Detector:  Since  tissue  phase-contrast  fringes  recorded  is  proportional  to  the 
Laplacian  and  gradient  of  phase-shifts,  hence  the  detector  should  have  high  sampling 
frequency.  We  employ  a  Konica  Minolta  Konica  24  x  30cm  photostimulable  phosphor  plate 
(CR-plate)  and  a  Konica  REGIUS  190  Reader  of  a  43.75pm  sampling  pitch. 

Image  Display  System:  For  the  high  resolution  mage-display  of  the  full  size  images  we 
will  use  the  GE  Seno  Advantage  review  workstation  (RWS)  in  the  Radiology  department  of 
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University  of  Alabama  at  Birmingham.  The  RWS  has  two  monitors  and  each  is  a  grayscale 
high-resolution  5  M  pixels  (2560  x2048)  CRT-monitor. 


(2).  Integrating  the  system  components  to  a  prototype  system 

Since  the  project  does  not  support  equipment  purchase,  most  of  the  system  components 
are  those  available  from  our  subcontractor  Dr.  Liu's  laboratory  at  University  of  Oklahoma.  Dr. 

Liu  group  performed  the  works  of  integrating  the  hardware  components  to  a  prototype  system. 

Mounting  the  System  on  An  Optical  Bench:  In  order  to  operate  the  prototype  in  both 
conventional  and  phase  contrast  imaging  modes  for  comparisons,  the  imaging  and 
measurement  components  (the  x-ray  tube,  CR-plate)  are  mounted  on  an  optical  rail,  which 
allows  to  vary  a  large  source-object  distance  R,  and  object-detector  distance  R2  as  needed.  It 
is  also  to  allow  use  of  a  precise  laser  alignment  process  for  later  measurements  of  x-ray  spectra 
and  system's  detective  quantum  efficiency  (DOE). 

The  source  is  a  micro-focus  tube  (Model  L8121-01 ,  Hamamatsu  Photonics)  with 
integrated  high  voltage  power  supply.  The  tube  has  a  tungsten-target  and  a  uniform  circular 
focal  spot  of  50|am  in  diameter.  The  tube  can  operate  at  150  kVp  with  maximal  tube  current  0.5 
mA.  During  the  hardware  integration,  the  tube  assembly  malfunction  occurred,  and  eventually 
the  tube  assembly  and  its  controller  failed.  After  tedious  efforts  in  dealing  with  the  OU  offices 
and  the  vendor  (Hamamatsu  Photonics)  for  the  repair  costs  and  repair  arrangements,  the 
source  assembly  is  currently  in  Hamamatsu's  shop  for  repair  and  calibrations.  As  a 
replacement  tube  for  the  time  being,  the  prototype  utilizes  an  micro-focus  tube  with  a 
Molybdenum  anode  and  a  Beryllium  output  window  (UltraBright  Microfocus  Source,  Oxford 
Instruments,  Scotts  Valley,  California).  The  x-ray  tube  can  be  operated  between  20  and  60  kVp 
and  between  10  and  60  W.  For  the  tests  performed  the  20  W-mode  was  utilized,  and  in  this 
mode  the  diameter  of  the  focal  spot  is  20  (mi .  Due  to  the  limitations  of  the  source,  the  high 
energy  x-ray  in  testing  was  that  from  the  tube  operating  at  60  kVp  only.  The  tube  voltage,  tube 
power,  and  exposure  time  are  controlled  precisely  by  a  software  application  interfacing  directly 
with  the  source. 

The  mounted  detector  is  a  Konica  Minolta  Konica  24  x  30cm  photostimulable  phosphor 
plate  (CR-plate)  and  a  Konica  REGIUS  190  Reader  of  a  43.75 (mi sampling  pitch. 

Prototype  Testing:  In  order  to  demonstrate  the  functionality  of  the  integrated  hardware 
component,  phantom  phase  contrast  imaging  experiment  were  performed. 

(a).  A  Tissue-Equivalent  Breast  Phantom:  The  Mammography  BR3D  phantom  (Model 
020,  CIRS,  Norfolk,  Virginia)  was  utilized  for  the  experiment.  The  BR3D  phantom,  as  a  tissue- 
equivalent  phantom,  was  fabricated  from  materials  simulating  100%  adipose  and  glandular 
tissues  blended  together  in  an  approximate  50/50  ratio  by  weight,  which  produces  a  tissue- 
equivalent  heterogeneous  background  on  an  x-ray  image.  The  phantom  consists  of  a  set  of  five 
pieces  of  the  same  size,  each  having  a  thickness  of  10  mm  and  a  different  blend  of  the 
materials.  The  combination  of  the  five  layers  was  designed  to  produce  realistic  tissue-equivalent 
x-ray  images.  The  qualitative  comparison  was  performed  for  a  high  energy  phase  contrast 
image  and  a  low  energy  conventional  radiograph.  The  ability  to  distinguish  structures  within  the 
images  will  provide  a  visual  indication  of  the  relative  image  quality.  The  conventional 
radiograph  of  the  tissue-equivalent  phantom  acquired  by  low  energy  (20  kVp)  conventional 


Page  9 


contact  mode  (M  =  1  ),  is  shown  in  Fig.  3-(left).  The  phase  contrast  image  of  the  phantom 
acquired  with  high  energy  (60  kVp)  and  R-\  =  R2  =  91.4  cm  is  shown  in  Fig.  3-(right).  The 
qualitative  comparison  illustrates  improved  tissue-material  contrast  in  the  phase  contrast  image 
as  compared  to  the  conventional  radiograph,  although  the  conventional  radiograph  was 
acquired  at  much  lower  tube  voltage  (20kVp)  and  it  would  have  more  attenuation-contrast 
compared  to  an  image  acquired  at  60  kVp.  The  testing  shows  the  prototype  functioning  well. 


Figure  3.  Acquired  Images  of  a  Mammography  BR3D  phantom,  (left):  the  image  acquired 
at  20  kVp,  20  W  in  conventional  contact  mode;  (right):  the  image  acquired  at  60  kVp,  20  W 
with  Rt  =  R2  =  91.4  cm  in  the  phase  contrast  mode. 


(b).  An  Acrylic  Edge  Phantom:  A  laser-cut  sharp  edge  of  a  thin  acrylic  plate  (~1 .5  mm) 
provides  a  phase  step:  a  phase  jump  from  acrylic  plate  to  air  along  the  sharp  edge.  The 
conventional  radiograph  of  the  edge  phantom  acquired  by  low  energy  (20  kVp)  conventional 
contact  mode  (M  =  1  ),  is  shown  in  Fig.  4-(left).  The  phase  contrast  image  of  the  phantom 
acquired  with  high  energy  (60  kVp)  and  in  phase  contrast  mode  is  shown  in  Fig.  4-(right). 
Compared  to  the  conventional  radiograph  of  the  acrylic  edge  shown  in  Fig.  4-(left),  the  phase 
contrast  image  of  the  same  phantom  at  60  kVp  (Fig.  4-(right))  shows  less  large-area  attenuation 
contrast  as  is  expected.  However,  the  phase  contrast  image  in  Fig.  4-(right)  uniquely 
demonstrates  the  signature  bright-dark  fringes  flanking  the  phase  step  edge.  The  diffraction 
fringes  is  proportional  to  the  Laplacian  of  the  phase  step.  In  fact,  the  intensity  profile  along  a 
orthogonal  line  across  the  edge  in  Fig.  4-(right)  demonstrates  a  downward  overshooting  and 
followed  by  an  upward  overshooting  in  the  phase  contrast  image  (Fig.  5-(right).).  But  such 
overshooting  disappears  in  the  intensity  profiles  of  conventional  radiograph  (Fig.  5-(left).) 
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Figure  4.  Acquired  Images  of  a  phase-step  edge  phantom,  (left):  the  image  acquired  at 
20  kVp,  20  W  in  conventional  contact  mode;  (right):  the  image  acquired  at  60  kVp,  20  W  in 
phase  contrast  mode  (  =  R2  =  91.4  cm.) 


Rows  Rows 


Figure  5.  (left):  the  intensity  profile  across  the  edge  for  the  conventional 
radiography  in  Fig.  4-(left);  (right):  the  intensity  profile  across  the  edge  for  the 
phase  contrast  image  in  Fig.  4-(right). 

In  summary,  the  system  hardware  integration  has  been  conducted  in  this  period.  A 
prototype  is  developed  and  the  tests  showed  that  the  prototype  is  functional  for  phase  contrast 
image  acquisition.  Once  the  in-repair  micro-focus  tube  (Model  L81 21  -01 ,  Hamamatsu 
Photonics)  is  returned  back  from  Hamamatsu  company,  we  will  soon  be  ready  to  acquire  phase 
contrast  images  at  120-150  kVp  as  designed.  In  addition,  we  will  continue  to  experiment  the 
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ways  to  match  the  CR-plate's  quantum  detection  efficiencies  with  the  designed  high  energy  x- 
rays  and  keep  reading  the  CR-plates  in  the  high  sampling  mode  of  the  Konica  REGIUS  190 
Reader. 


Key  Research  Accomplishments 


•  Phase  retrieval  from  a  single  phase  contrast  image  is  the  key  strategy  for  this  project.  In 
this  period  we  successfully  developed  a  new  phase  retrieval  method  with  a  full  account 
of  the  nonlinear  Fresnel  propagators.  This  new  phase  retrieval  method  provides  the 
scientific  foundation  and  framework  for  developing  new  and  more  accurate  phase 
retrieval  algorithms  for  this  project.  The  computer  simulation  demonstrated  superiority  of 
this  new  method  over  our  previous  phase-retrieval  method  in  reconstructing  the  high 
frequency  components  of  imaged  objects. 

•  Conducted  the  system  design  study  for  the  project.  Especially  we  have  determined 
the  optimal  imaging  geometry  based  on  the  quantitative  phase  contrast  visibility 
analysis. 

•  Conducted  system  hardware  integration.  A  prototype  is  developed  and  the  tests 
showed  that  the  prototype  is  well  functioning  for  phase  contrast  image  acquisition. 


Reportable  Outcomes 


In  this  project  period  we  have  published  the  research  results  in  peer-reviewed  journals 
and  presented  the  work  in  an  international  conferences,  as  listed  in  the  following. 

Peer-Reviewed  Journal  Article: 

X.  Wu  and  A.  Yan:  "Phase  retrieval  from  one  single  phase  contrast  x-ray  image,"  Optics 
Express  17:  11187-11196  (2009). 

This  journal  is  the  top-ranking  journal  in  terms  of  the  Impact  Factor  among  all  sixty-four 
journals  in  the  field  of  Optics. 


Published  Abstract  and  Conference  Presentation: 

X.  Wu,  A.  Yan  ,  H.  Liu:  “Contrast  mechanism  and  origin  of  artifacts  in  phase-sensitive 
x-ray  volumetric  imaging”,  Medical  Physics  36:  2785  (2009). 

X.  Wu,  A.  Yan  ,  H.  Liu:  “Contrast  mechanism  and  origin  of  artifacts  in  phase-sensitive  x- 
ray  volumetric  imaging.”  Oral  Presentation  at  the  51st  Annual  Meeting  of  the  American 
Association  of  Physicists  in  Medicine,  July  29,  2009,  Anaheim,  California 
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Conclusion 


With  the  support  by  USAMRMC  we  have  successfully  conducted  studies  on  x-ray  phase 
imaging  in  the  first  period  of  this  project.  As  the  key  strategy  of  the  project  is  to  perform  the 
phase  retrieval  from  a  single  phase  contrast  image,  in  this  period  we  successfully  developed  a 
new  phase  retrieval  method  from  a  single  image  with  a  full  account  of  the  nonlinear  Fresnel 
propagators.  This  new  phase  retrieval  method  provides  the  scientific  foundation  and 
framework  for  developing  new  and  more  accurate  phase  retrieval  algorithms  for  this  project. 

The  computer  simulation  demonstrated  superiority  of  this  new  method  over  our  previous  phase- 
retrieval  method  in  reconstructing  the  high  frequency  components  of  imaged  objects. 

In  this  period  we  have  conducted  the  system  design  study.  The  crucial  focus  of  the 
system  design  is  to  determine  the  imaging-geometry,  which  should  optimize  the  phase  contrast 
visibility  by  balancing  the  conflict  requirements  of  x-ray  spatial  coherence,  large  diffraction 
fringes  while  being  constrained  by  the  finite  sizes  of  the  tube  focal  spot  and  detector-pitch.  We 
investigated  how  the  relative  phase-visibility  factor  varies  with  system  hardware  parameters  and 
found  the  optimal  reference  imaging  geometry  parameters  for  the  project.  To  implement  the 
system  design,  a  prototype  has  been  built  by  integrating  the  system  components.  With  the 
system  prototype  phase  contrast  images  were  acquired  with  different  phantoms  for  verifying  the 
prototype's  functionality. 

In  the  coming  year  we  will  continue  the  efforts  in  improving  the  prototype's  performance, 
especially  work  with  the  repaired  tube  operating  on  the  designed  high  tube  voltages,  and  with 
the  matched  CR  plates  readable  with  the  high  sampling  mode.  We  will  conduct  phase  imaging 
experiments  with  custom-made  breast  phantoms,  and  will  reconstruct  the  phase  maps  of  the 
phantoms.  Especially  we  will  incorporate  the  measured  source  incoherence  effects  and  the 
detector  response  effects  into  the  phase  retrieval  algorithms  for  accurate  phase  retrieval.  We 
will  start  to  experimentally  verify  the  accuracies  of  the  reconstructed  maps  of  the  projected 
electron  densities  of  the  imaged  phantoms. 


References 


Wu,  X.,  Liu,  H.  2004.  A  new  theory  of  phase-contrast  x-ray  imaging  based  on  Wigner 
distributions.  Medical  Physics  31 , 2378-2384. 

Wu,  X.,  Liu,  H.,  Yan,  A.  2005.  X-Ray  phase-attenuation  duality  and  phase  retrieval. 
Optics  Letters  30,  379-381 . 

Wu,  X.,  Yan,  A.  2009.  Phase  retrieval  from  one  single  phase  contrast  x-ray  image, 
Optics  Express  17,  11187-11196. 


Page  13 


Appendices 


Copies  of  following  publications  are  enclosed  for  review. 
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x-ray  volumetric  imaging”,  Medical  Physics  36:  2785  (2009). 
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Abstract:  Phase  retrieval  is  required  for  achieving  artifact-free  x-ray 

phase-sensitive  3D  imaging.  A  phase-retrieval  approach  based  on  the 
phase-attenuation  duality  with  high  energy  x-rays  can  greatly  facilitate 
for  phase  sensitive  imaging  by  allowing  phase  retrieval  from  only  one 
single  projection  image.  The  previously  derived  phase  retrieval  formula  is 
valid  only  for  small  Fresnel  propagator  phases  corresponding  to  common 
clinical  imaging  tasks.  In  this  work  we  presented  a  new  duality-based  phase 
retrieval  formula  that  can  be  applied  for  cases  with  large  Fresnel  propagator 
phases  corresponding  to  high  spatial  resolution  imaging.  The  computer 
simulation  demonstrated  superiority  of  this  new  formula  over  the  previous 
phase  retrieval  formula  in  reconstructing  the  high  frequency  components 
of  imaged  objects.  A  modified  Tikhonov  regularization  technique  has 
been  devised  for  phase  retrieval  in  cases  of  very  high  resolution  and  large 
object-detector  distance  such  that  some  Fresnel  propagator  phases  may  be 
close  or  greater  than  n.  This  new  phase  retrieval  formula  lays  the  foundation 
for  implementing  high-resolution  phase-sensitive  3D  imaging  of  soft  tissue 
objects. 

©  2009  Optical  Society  of  America 
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1.  Introduction 

Conventional  x-ray  imaging  for  soft  tissues  (such  as  breast,  brain,  liver,  etc.)  is  limited  in  its  sen¬ 
sitivity  for  detecting  subtle  tissues  pathological  changes,  since  the  imaging  relies  on  the  small 
difference  in  x-ray  attenuation  between  the  lesions  and  soft  tissues  of  variable  structure.  How¬ 
ever,  as  x-ray  wave  passes  through  tissues,  x-ray  undergoes  phase-shifts  as  well.  The  amount 
of  x-ray  phase  shift  (j)  by  tissues  is 


where  re  is  the  classical  electron  radius,  h  is  the  Plank  constant  and  c  the  speed  of  light,  E  is 
x-ray  photon  energy.  Here  pe  denotes  tissue  electron  density,  and  the  integration  of  electron 
density  pe  is  over  the  x-ray  path  and  the  integral  is  the  projected  electron  density  pe  p  [1,  2,  3]. 
It  turns  out  that  x-ray  phase-shift  differences  between  tissue  and  lesions  are  about  one  thou¬ 
sand  times  larger  than  their  attenuation  differences  [1,2,  3].  Hence  the  phase-sensitive  imaging 
has  the  potential  to  greatly  enhance  the  lesion  detection  sensitivity  and  specificity,  and  reduce 
the  radiation  dose  associated  with  x-ray  imaging.  In  the  inline  phase  contrast  imaging  as  stud¬ 
ied  in  previous  works,  the  x-rays  with  shifted  phases  diffract  from  the  object  to  the  detector, 
and  the  diffracted  x-rays  form  bright  and  dark  fringes  at  tissue  boundaries.  The  edge  enhance¬ 
ment  thus  generated  relies  on  the  spatial  coherence  of  the  x-ray  source,  and  Laplacian  and 
gradients  of  x-ray  phase-shifts  caused  by  the  object,  and  gradients  of  the  objects  attenuation 
[4,  5,  6,  7,  8,  9,  10,  11,  12,  13,  14].  The  procedure  of  disentangling  tissue  phase-shifts  from 
the  mixed  contrast  mechanism  and  retrieving  the  phase  maps  from  acquired  phase  contrast 
images  is  called  the  phase  retrieval.  Phase  retrieval  technique  plays  a  central  role  in  phase- 
sensitive  x-ray  imaging.  By  means  of  phase  retrieval,  one  can  reconstruct  a  quantitative  map 
of  phase-shifts,  a  phase  image  of  the  imaged  object  [3,  6,  13,  15,  16].  According  to  Eq.  (1), 
a  retrieved  phase  image  is  equivalently  a  map  of  imaged  objects  quantitative  projected  elec¬ 
tron  densities.  Therefore,  tissue  phase  images  can  provide  quantitative  tissue  characterization. 
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which  potentially  can  increase  the  sensitivity  and  specificity  of  a  diagnosis.  Moreover,  phase 
retrieval  is  a  necessary  procedure  for  phase-sensitive  volumetric  imaging  such  as  the  phase- 
sensitive  tomography  and  tomosynthesis  [7,  14,  15].  In  the  phase-sensitive  imaging,  x-ray  does 
NOT  propagate  along  straight  lines  because  of  x-ray  refraction  and  diffraction  at  tissue  inter¬ 
face  and  boundaries.  The  diffraction  fringes  generated  by  phase  contrast  at  tissue  interfaces 
will  obviously  render  conventional  volumetric  image  reconstruction  methods  erroneous,  since 
conventional  volumetric  reconstruction  methods  assume  that  x-rays  propagate  along  straight 
lines.  The  phase-sensitive  CT  experiments  showed  that  conventional  tomography  reconstruc¬ 
tion  algorithms  applied  to  phase-contrast  projections  can  lead  to  erroneous  tomography  images 
with  artifacts  such  as  the  streaks  and  negative  densities  [7].  The  exact  mechanism  of  the  phe¬ 
nomenon  was  revealed  in  one  of  our  recent  work.  In  fact,  if  we  denote  the  imaged  objects 
linear  attenuation  coefficient  by  p(r)  and  its  electron  density  by  pe(r),  with  the  conventional 
filtered  backprojection  method  the  reconstructed  “apparent  attenuation  coefficients”  /ireCon (/) 
was  found  as  [17]: 


(2) 


where  R2  denotes  the  object-detector  distance,  A  the  x-ray  wavelength.  Eq.  (2)  shows  that  the 
reconstructed  tomogram  consists  of  three  sub-tomograms.  The  first  sub-tomogram  is  the  tomo¬ 
gram  of  the  tissue  linear  attenuation  coefficients  p(r).  The  second  sub-tomogram  is  a  scaled 
map  of  the  3-D  Laplacian  of  tissue  electron  densities.  The  third  sub-tomogram  /jmiXed(Vj  repre¬ 
sents  the  artifacts  reflecting  the  global  variation  of  both  /r(r)  and  p,, ( r) ,  but  is  not  representing 
the  local  variations  of  ju(r)  and  pr.  ( r) .  It  is  this  term  that  causes  artifacts  such  as  the  streaks  in 
reconstruction  images.  In  order  to  correctly  deal  with  volumetric  reconstruction  issue,  as  a  gen¬ 
eral  strategy  one  should  first  disentangle  tissue  phase-shifts  from  the  mixed  contrast  mechanism 
from  acquired  phase-sensitive  projections  and  retrieve  the  phase  maps  for  all  projections,  then 
reconstruct  the  3D  images  of  object  attenuation  coefficients  and  electron  densities,  respectively. 

For  phase  retrieval  one  in  general  needs  to  acquire  at  least  two  and  mostly  four  to  five  im¬ 
ages  with  different  object-detector  distances  [3,  15].  This  requirement  is  imposed  by  the  x-ray 
wave  nature  and  its  wave  equation  [3,  8,  9,  10].  Multiple  image-acquisitions  for  a  phase  re¬ 
trieval  increase  radiation  dose  associated  with  the  technique.  For  phase  sensitive  tomography  or 
tomosynthesis,  multiple  image-acquisitions  for  each  projection-angle  make  phase-sensitive  to¬ 
mography  and  tomosynthesis  implementation  especially  cumbersome:  multiple  phase-contrast 
images  have  to  be  acquired  with  different  object-detector  distances  for  all  the  projections  dur¬ 
ing  a  scan,  and  the  number  of  required  projections  can  be  as  high  as  one  thousand.  In  addition, 
the  image  data  alignment  procedure  is  quite  time-consuming.  Moreover,  radiation  dose  is  mul¬ 
tiplied  for  multiple  image  acquisition  per  projection  view  as  well.  Therefore,  an  important  chal¬ 
lenge  is  to  find  an  effective  and  low-dose  approach  for  phase  retrievals  for  phase-sensitive  imag¬ 
ing.  Many  studies  were  conducted  addressing  this  challenge  [6,  13,  14].  We  proposed  a  unique 
solution  for  phase-sensitive  imaging  with  high  energy  x-rays.  We  noted  that  with  increasing  x- 
ray  energy  tissues  photoelectric  absorption  decreases  as  1  /En  with  rk,  3—4.  but  tissues  phase 
shift  decreases  much  slower,  and  only  as  l/E  as  is  shown  in  Eq.  (1).  Therefore,  phase-sensitive 
imaging  with  x-rays  can  fully  take  advantage  of  phase  contrast  and  greatly  reduce  radiation 
doses  associated  with  x-ray  imaging,  by  taking  advantage  of  the  dual  relation  between  the  high 
energy  x-rays  phase-shifts  and  attenuation  for  soft  tissues.  Furthermore,  analyzing  the  tissue 
attenuation  data,  we  observed  that  the  soft  tissue  attenuation  cross  sections  are  very  well  ap¬ 
proximate  by  that  of  x-ray  incoherent  scattering  for  x-rays  of  about  60  keV  to  500  keV.  Under 
this  circumstance,  we  noted  that  soft-tissue  attenuation  and  soft-tissue  phase  are  all  related  to 
the  projected  electron  density,  we  called  this  new  notion  as  the  phase-attenuation  duality.  Based 
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Fig.  1.  A  shematic  diagram  for  the  measuring  system 


on  this  concept  of  phase-attenuation  duality,  we  found  a  way  to  retrieve  a  phase-map  from  only 
one  single  recorded  image,  that  is,  the  retrieved  phase  map  is  given  by  [13]: 


(3) 


where  Ri  is  the  distance  from  the  point  x-ray  source  to  object,  AS  is  the  distance  from  object 
to  detector,  M  =  (R\  -I  Ri) / R\  is  the  geometric  magnification,  r  is  the  position  vector  (see 
Fig.  1.)  In  Eq.  3,  I(Mr\R\  +  Rt)  is  the  measured  image  intensity  at  detector  plane,  and  Im  is 
the  intensity  at  entrance,  the  operator  V2  denotes  the  Laplacian  operator.  In  Eq.  (3)  Okn  is  the 
Klein-Nishina  total  cross  section,  that  is,  the  total  cross  section  for  x-ray  Compton  scattering 
with  a  free  electron  [17].  Note  that  Okn  is  x-ray  energy  dependent,  as  is  shown  below.  We  point 
out  that  the  high  sensitivity  of  x-ray  phase  change  is  realized  as  the  high  ratio  Xrej 1 
in  this  circumstance.  In  fact,  for  clinical  imaging  photon  energy  is  no  higher  than  150  keV, 
and  in  this  case  Xrej Okn  ~  103.  In  other  words,  with  this  duality  the  high  sensitivity  of  phase 
imaging  manifests  itself  as  the  high  ratio  of  x-ray  wavelength  over  the  classical  electron  radius. 
Compared  to  the  common  multiple-acquisition  based  phase  retrieval  approaches  in  the  literature 
[3,  6,  15],  our  duality-based  phase  retrieval  method  is  a  conveniently  implemented  method 
with  tremendous  dose-saving  for  many  potential  applications  such  as  mammography.  These 
advantages  will  be  critical  for  phase-sensitive  volumetric  imaging  such  as  breast  tomosynthesis 
and  computed  tomography. 

However,  there  is  a  limitation  on  applicability  of  the  phase  retrieval  formula  Eq.  (3).  In  fact. 
Eq.  (3)  is  not  valid  for  the  cases  with  very  high  spatial  resolution  and  with  large  object-detector 
distances.  The  physics  behind  this  can  be  intuitively  explained  as  follows.  In  fact,  in  inline 
phase-sensitive  imaging  the  phase  contrast  is  formed  by  the  diffraction  from  the  object  to  the 
detector  of  x-rays  with  shifted  phases.  The  Fresnel  diffraction  can  be  thought  as  a  x-ray  wave- 
front  filtering  with  the  Fresnel  propagator  exp(inXzu2),  where  X  is  the  wavelength  and  u  is 
the  wave-front  transverse  spatial  frequency  vector,  and  z  =  Rj/M  in  terms  of  the  parameters  in 
Eq.  (3).  So  the  phase  of  the  propagator,  i.e.  nXzu2,  determines  how  strongly  the  wave-front  exit 
from  the  object  will  be  diffracted  in  x-ray  propagation.  The  phase  retrieval  formula  Eq.  (3)  was 
derived  by  linearizing  Fresnel  propagator  phase  under  the  assumption  of  nXzu 2  <C  1  [8,  13]. 
This  linearizing  is  equivalent  to  simplifying  x-ray  wave  equation  of  motion  to  the  transport  of 
intensity  of  equation  (TIE)  [8].  However,  for  cases  with  very  high  spatial  resolution  and  with 


#1 10095  -  $15.00  USD  Received  14  Apr  2009;  revised  20  May  2009;  accepted  17  Jun  2009;  published  19  Jun  2009 


large  object-detector  distances,  such  as  cases  with  phase-sensitive  imaging  experiments  with 
synchrotron  based  x-ray  sources,  or  the  high  resolution  imaging  experiments  with  biological 
specimen  or  small  animals,  the  high  spatial  resolution  or  large  propagation  distances  make 
Fresnel  propagator  phases  nXzu1  larger  than  1  for  many  transverse  spatial  frequencies  of  x-ray 
wave-front.  For  example,  for  60  keV  x-rays,  if  M  =  1  and  AS  =  2m,  detector  pixel  is  5  /im,  then 
nXzu 2  =  1.299.  For  these  cases  one  cannot  linearize  Fresnel  propagators  and  should  consider 
the  full  exponential  propagators  instead.  In  section  2,  we  will  first  briefly  review  the  concept 
of  the  phase-attenuation  duality  and  how  it  greatly  simplify  the  phase  retrieval  for  soft-tissue 
objects.  We  will  then  derive  a  new  and  general  PA-Duality  based  phase  retrieval  formula,  as  the 
main  result  of  this  work.  With  this  general  formula,  one  can  retrieve  the  phase  from  one  single 
phase  contrast  image  and  is  valid  for  cases  involved  with  large  Fresnel  propagator  phases. 
Hence  this  new  phase  retrieval  formula  extends  its  applicability  to  cases  with  very  high  spatial 
resolution  and  with  large  object-detector  distances.  In  section  3,  we  show  the  simulation  results 
with  the  general  formula  for  phase  retrieval. 


2.  Phase-attenuation  duality  and  phase  retrieval  from  one  single  phase  contrast  image 
for  soft  tissue  objects 

As  is  explained  in  previous  works  [3,  4,  8,  9,  10],  in  phase-sensitive  imaging  the  imaged  object 
can  be  modeled  by  a  two-dimensional  transmission  function  T{r)  =  Ao(r)el^r\  where  0o(r) 
denotes  the  x-ray  phase-shift  caused  by  the  object,  and  Ao(?)  denotes  the  x-ray  transmission, 
or  the  attenuation-map  of  object.  For  sake  of  clarity  in  exploring  the  effects  of  non-linearity  of 
the  Fresnel  propagators,  we  assume  the  x-ray  source  is  a  quasi-monochromatic  point  source  in 
following  derivations.  We  will  discuss  more  general  x-ray  sources  in  section  3. 

Let  R\  denote  the  source-object  distance  and  AS  the  object-detector  distance.  The 
geometric  magnification  factor  M  is  equal  to  (  R\  +  AS) / R  \  .  We  also  assume  that 
one  encounters  only  moderate  variation  of  the  object  phase  function  0(r)  such  that 
|0  (r  +  XRzu/lM)  —  0  {f  —  XRju j2M)\  <C  1  [3].  Transmitting  from  the  object,  x-rays  undergo 
attenuation  and  phase-shift,  and  diffract  over  a  distance  R2  to  the  detector.  An  important  task  is 
to  find  out  the  formula  of  the  image  intensity  as  function  of  the  objects  transmission  function 
and  imaging  geometry.  A  general  solution  had  been  found  in  one  of  our  previous  works  [3]. 
Starting  from  the  paraxial  Fresnel-Kirchhoff  diffraction  theory,  we  found  out  in  this  work  that 
the  Fourier  transform  of  x-ray  irradiance  at  detector  plane  is  given  by: 
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where  I(jd)  denotes  the  x-ray  irradiance  at  the  detector  entrance,  and  the  symbol  -IX  (■)  denotes 
the  2-D  Fourier  transform  and  Im  is  the  intensity  of  the  incident  x-ray  upon  the  object,  and 
u  is  the  spatial  frequency  vector  in  the  object  plane.  Note  that  the  Fresnel  propagator  phases 
nXRiu2 /M  enter  Eq.  (4)  as  the  argument  of  sinusoidal  transfer  functions.  For  cases  with  mod¬ 
erate  spatial  resolutions  u  and  AS  ,  Eq.  (4)  is  much  simplified,  and  the  inverse  Fourier  transfor¬ 
mation  of  the  simplified  Eq.  (4)  can  be  written  as  [8]: 
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Equation  (5)  is  the  TIE-based  formula  for  image  irradiance,  which  was  first  derived  by  Nugent 
and  colleagues [3].  Recently  Guigay  et  al  showed  the  limitations  of  Eq.  (5)  with  theoretical 
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analysis  and  experiments  [17].  They  found  that  for  in-line  phase-sensitive  imaging  with  long 
R2  and  high  resolutions  such  that  nXR2u2  /M  >  1  ,  the  phase  retrieval  based  on  Eq.  (5)  would  be 
inaccurate,  one  has  to  use  the  general  formula  Eq.  (4)  for  accurate  quantitative  phase  imaging. 

Furthermore,  we  note  that  soft  tissues  encountered  in  clinical  imaging  are  composed  almost 
exclusively  by  atoms  of  the  light  elements  with  the  atomic  numbers  Z  <  10.  For  example,  the 
composition  of  breast  tissue  consists  mainly  of  light  elements  such  as  hydrogen,  carbon,  nitro¬ 
gen  and  oxygen.  In  contrast,  the  sum  of  weight  fractions  of  other  heavier  elements  in  breast  tis¬ 
sue  is  only  less  than  1%  [18].  The  same  is  true  for  other  soft  tissues  such  as  the  brain  grey/white 
matter  [18].  Analyzing  the  tissue  attenuation  data,  we  observed  that  the  soft  tissue  attenuation 
cross  sections  are  very  well  approximated  by  that  of  x-ray  incoherent  scattering  for  x-rays  of 
about  60keV  to  500keV[19].  In  addition,  the  incoherent  scattering  functions  become  linearized 
as  well  [19].  Under  this  circumstance,  while  the  phase-shift  is  caused  by  the  coherent  x-ray 
scattering,  but  soft  tissue  attenuation  is  determined  by  the  incoherent  x-ray  scattering  for  soft 
tissues,  so  both  the  x-ray  phase-shift  and  attenuation  by  soft  tissues  are  all  determined  by  the 
projected  electron  density  for  these  high  energy  x-rays.  We  called  this  complementary  relation¬ 
ship  between  phase-shift  and  attenuation  for  soft  tissues  as  the  phase-attenuation  duality.  When 
the  P-A  duality  holds,  x-ray  attenuation  and  phase  shift  by  the  object  is  related  in  following 
way: 

Hr)  =  -ArfpejP(r),  A20(r)  =  exp(-OkNpe,„(r)),  (6) 

where  Okn  is  the  total  cross  section  for  x-ray  Compton  scattering  with  a  free  electron: 
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where  rj  =  Ephoton/mec2,  here  we  denote  the  photon  energy  of  the  primary  x-ray  beam  by 
/-photon,  and  mec2  is  the  rest  electron  energy  and  equal  to  51  IkeV.  The  duality  condition  can  be 
used  to  simplify  Eq.  (4).  In  fact,  we  can  rewrite  Eq.  (4)  as 
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Especially  with  the  duality  condition  we  found 
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In  Eq.  (9)  we  made  an  approximation  of  dropping  the  negligible  contribution  of  the  higher 
order  terms  of  Okn/^tv  The  approximation  is  based  on  GknMt>  ~  10  and  we  assumed 
|p  (r  +  XR2u/2M)  —  0  (r  —  XR2u/2M)\  -C  1,  moderate  variations  of  object  phase  shifts.  Sub¬ 
stituting  Eqs.  (9)  into  Eq.  (8),  we  found: 
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Equation  (10)  is  the  central  result  of  this  work,  it  is  the  formula  for  x-ray  irradiance  for  high 
resolution  soft  tissue  imaging  with  high  energy  x-rays.  In  the  derivation  the  parameter  A re! Okn 
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plays  a  significant  role.  It  is  interesting  to  note  that  Gureyev  et  al  obtained  a  similar  result  in 
ref.[20]  (Eq.(24))  with  the  exception  that  the  term  nXRiu1  /  M  in  Eq.  (10)  was  absent.  The  dif¬ 
ferences  may  lie  in  the  different  assumptions  made  regarding  the  object  transmission  functions 
and  x-ray  energies.  Gureyev  and  coauthors  assumed  the  objects  as  single  material  homoge¬ 
neous  objects  [20],  but  our  results  are  for  heterogeneous  objects  of  the  light  elements  imaged 
with  high  energy  x-rays  (60  keV  or  higher),  i.e.,  for  the  cases  the  phase-attenuation  duality 
holds.  In  addition,  in  their  derivation  they  separated  the  slowly  and  rapidly  varying  components 
of  the  amplitude  and  phase,  and  they  assumed  that  the  magnitude  of  the  rapidly  varying  com¬ 
ponents  are  much  smaller  than  1,  but  we  did  not  make  such  limiting  assumption  in  derivation  of 
Eq.  (10).  We  stress  that  the  Klein-Nishina  cross-section  is  energy  dependent  shown  in  Eq.  (7) 
and  X  re  j (Jkn  is  about  1000  for  inline  x-ray  energy  range  used  in  diagnostic  imaging.  This  large 
Are/<TKN  represents  the  high  sensitivity  of  phase  contrast  in  tissue  differentiation,  but  it  also 
makes  the  x-ray  irradiance  formula  Eq.  (10)  getting  very  much  simplified.  It  is  easy  to  verify 
that  for  the  low  resolution  cases  such  as  those  with  n'kR^u1  jM  <C  1,  Eq.  (10)  is  reduced  to 


we  recover  our  previous  result  as  shown  in  [13].  On  the  other  hand,  for  general  cases 
nXR^u2  /M  ~  1  or  nXRiu1  /M  >  1  Eq.  (11)  cannot  be  applied,  and  one  should  use  Eq.  (10) 
for  phase  retrieval  as  is  shown  below. 

3.  Simulation  results 

Comparing  Eq.  (10)  and  Eq.  (11)  for  use  in  performing  phase  retrieval,  the  only  difference  lies 
in  the  inverse  filters  derived  for  retrieving  JF  (Aq).  For  phase  retrieval  based  on  Eq.  (10)  the 
inverse  filter  is  found  as 
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where  a  =  tiXRju1  j  M ,  the  FT-space  Fresnel  Propagator  phase  as  defined  in  section  1,  and 
y  =  X  re  / <7k x •  which  represents  the  phase  contrast  enhancement  as  explained  earlier.  8  = 
1  T  (2y+  a)2  and  ai  =  a  —  sin-1  [(2y+  a) /8\.  On  the  other  hand,  for  phase  retrieval  based 
on  Eq.  (11)  the  inverse  filter  will  be 


1 

l  +  2ya' 


(13) 


Since  y  ~  103,  the  filters  in  Eq.  (12)  and  Eq.  (13)  are  close  when  a  is  small.  When  a  is  getting 
bigger,  close  to  or  greater  than  n  (oq  close  to  n/2),  the  difference  between  the  two  inverse 
filters  is  getting  bigger.  Figure  2  shows  a  plot  of  the  two  filters,  with  respect  to  the  change  of  a. 
The  inverse  filter  in  Eq.  (12),  depicted  in  Fig.  2  as  the  solid  line,  changes  the  sign  when  a  passes 
through  kn  (and  thus  a.\  passes  through  (k—  \/2)n)  and  diverges  to  °°  when  ai  =  (k—  1/2 )n, 
k=  1,2,  •  •  • .  On  the  other  hand,  the  filter  in  Eq.  (13),  depicted  as  the  dashed  line  in  Fig.  2,  keeps 
positive  and  approaches  zero  as  a  gets  bigger.  So  the  difference  of  the  retrieved  phase  maps 
using  Eq.  (10)  and  Eq.  (11)  will  be  highly  dependent  on  the  value  of  the  Fresnel  propagator 
phase  a  =  nXR2u2/M. 

To  demonstrate  this  phenomenon,  we  performed  two  simulations.  First,  we  set  the  phase  shift 
(j>  to  be  a  sine-function  in  x  direction  with  period  of  32  pixel  size  and  constant  in  y  direction. 
Assuming  the  pixel  pitch  is  p  jum,  the  spatial  frequency  corresponding  to  0,  in  x  direction,  is 
RrO  =  l/32y>  1/jum.  In  our  simulations,  for  the  purpose  of  simplification,  the  magnification  M 
as  well  as  the  intensity  at  entrance,  /m,  are  assumed  to  be  1.  The  x-ray  photon  energy  is  set 
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Fig.  2.  Plots  of  filters  (13)  and  (12) 


to  60  keV.  This  gives  the  x-ray  wavelength  A  =  2.0667  x  10  11  m  and  the  total  cross  section 
for  x-ray  Compton  scattering  with  a  free  electron  <7kn  =  5.4526  x  1 0  29  m2  as  determined 
from  Eq.  (7).  The  classical  electron  radius  re  is  set  to  2.82  x  10~15  m.  By  assuming  the  de¬ 
tectors  resolution  p  =  1  /im  and  the  phase  shift  <l>(x.y)  =  0.5sin(32/?x)  —  1,  according  to  the 
phase-attenuation  duality  discussed  earlier,  the  attenuation  map  is  related  to  the  phase  map 
as  A/(x,y)  =  exp  [OYm/Xre  •</>]  =  exp  [(sin(32px)  —  2) /2y],  with  y=  Aty/okN  =  1-0726  x  103. 
We  then  performed  Fresnel  diffraction  to  determine  the  phase-contrast  image  irradiances  at  dis¬ 
tances  downstream [21],  and  the  retrieved  maps  of  4/  were  found  by  using  the  inverse  filters  of 
Eq.  (12)  and  Eq.  (13),  which  were  derived  based  on  Eq.  (10)  and  Eq.  (11),  respectively.  Because 
of  the  large  value  of  y,  phase  contrast  image  irradiances  {/}  were  found  be  maps  with  each  I  be¬ 
ing  dominated  by  two  single  frequency  component  at  ux  =  ±;/aq  (uxq  =  1  /32 p  =  1/32  x  106  m) 
plus  a  zero  frequency  component.  This  being  so,  the  shape  and  variation  of  the  retrieved  Ag- 
maps  were  determined  by  their  dominant  frequency  component  at  frequencies  ux  =  ±ux 0.  In  or¬ 
der  to  clearly  show  the  differences  between  phase  retrievals  based  on  Eq.  (10)  and  Eq.  (1 1),  we 
examined  the  ratio  of  the  dominant  frequency  components  of  the  two  retrieved  Ajj-maps  based 
on  Eq.  (10)  and  Eq.  (11),  respectively.  We  found  the  ratio  being  r  =  (1  +2ya)/(<5  -cos(ai)) 
for  a  given  object-detector  distance  Ri-  This  ratio  reflects  the  capability  of  achieving  accu¬ 
rate  phase  retrieval  with  Eq.  (10)  against  Eq.  (11).  For  small  1,  that  is,  for  small  Fres¬ 
nel  propagator  phases,  this  ratio  is  1,  so  Eq.  (10)  against  Eq.  (11)  both  are  equally  good  for 
the  accurate  phase  retrieval.  However,  when  a  is  chosen  such  that  cosjcq)  is  close  to  0,  this 
ratio  can  be  very  large.  For  example  by  setting  R2  =  49.55  m,  we  have  a  mentioned  above 
equals  to  n,  the  ratio  r  =  (1  +2y)/(t>  ■  cos(ai))  ss  6700.  This  large  ratio  means  that  the  re¬ 
trieved  ,4 /-map  based  on  Eq.  (11)  has  diminished  image-contrast  and  erronous  values,  as  is 
shown  in  Fig. 3(c).  The  phase  retrieval  result  based  on  Eq.  (10),  shown  in  Fig.  3(d),  how¬ 
ever  is  much  more  accurate.  When  R2  =  49.55  m,  the  FT-space  Fresnel  Propagator  phase 
a  =  nXRyu2 /M  ss  n.  So  we  have,  from  Eq.  (10),  $  ((M2//in)/)  ss  —  &  (Aq),  which  is  differ¬ 
ent  from  the  result  &  ((M2 //;„)/)  «  [1  +27r(Are/oKN)]  &  (Ag)  derived  from  the  conventional 
TIE  prediction  (Eq.  (11)),  which  does  not  hold  for  such  large  Fresnel  Propagator  phases.  This 
phenomenon  is  shown  in  Fig.  (3)(b),  where  the  intensity  I  is  acquired  with  Fresnel  diffraction, 
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in  which  the  dark-white  area  in  A5  is  reversed.  Note  that  the  inverse  filter  in  Eq.  (12)  is  singular 
at  oq  =  [k—  1/2)  tc,  k  =  1 ,2,  •  •  • ,  because  of  the  zero-crossing  of  cos(ai).  To  deal  with  the  sin¬ 
gularity  associated  with  the  zero-crossing  of  cos(ai),  one  should  regularize  the  inverse  filter. 
A  common  regularization  technique  is  the  Tikhonov  regularization,  which  modifies  1  / cos(ai) 
in  the  filter  of  Eq.  (12)  to  cos(ai)/(cos(ai)2  +  jc2),  where  K  is  the  Tikhonov  regularization 
parameter.  In  this  specific  example  of  the  sinusoidal  object,  however,  since  we  know  &  (/)  is 
dominated  by  frequency-components  only  at  ux  =  ±ux o  and  ( u ,  v)  =  (0, 0),  by  setting  &  (/)  =  0 
at  all  other  frequency  points  purposely,  we  can  avoid  the  problem  caused  by  the  zero-crossing 
of  the  inverse  filter. 


(a)  (b)  (c)  (d) 

Fig.  3.  Comparation  of  the  quality  of  the  retrieved  results  Aq  using  Eq.  (11)  and  Eq.  (10). 
A5  =  exp[0/y],  where  phase  shift  (j)  =  0.5sin(32px)  —  1,  pixel  pitch  p=  1pm,  and  object  to 
detector  distance  R2  =  49.55m.  (a)  is  the  true  Ajp  (b)  is  the  phase-contrast  intensity  image 
I  acquired  from  Fresnel  diffraction;  (c)  is  retrieved  A(2  using  Eq.  (11);  (d)  is  retrieved 
using  Eq.  (10). 


In  another  simulation,  we  set  all  the  parameters  the  same  as  the  previous  simulation  except 
replacing  the  test-object  with  the  well-known  Lena  image  as  is  shown  in  Fig.  4(a).  The  simula¬ 
tion  results  for  the  case  of  R2  =  5  m  are  shown  in  Figures  4(c)-4(d).  The  phase-contrast  image 
I,  shown  in  Fig.  4(b),  is  acquired  with  Fresnel  diffraction  from  the  object  downstream  to  the 
detector  [21].  Figure  4(c)  shows  the  retrieved  image  with  the  inverse  filter  of  Eq.  (13),  which 
is  based  on  Eq.  (11).  Since  the  filter  of  Eq.  (13)  is  incorrect  when  the  Fresnel  propagator  phase 
a  is  big,  say  close  or  greater  than  n,  the  high  frequency  components  of  the  object  was  unduly 
suppressed  in  the  phase-retrieval.  This  deficiency  is  reflected  in  Fig.  4(c)  such  that  the  retrieved 
image  gets  blurred  in  locations  with  high  contrast  changes,  such  as  the  hair,  the  edge  of  the  hat 
and  body.  In  contrast,  the  phase  retrieval  using  the  inverse  filter  of  Eq.  (12),  which  is  based  on 
Eq.  (10),  preserves  objects  high  frequency  components,  as  is  shown  in  Fig.  4(d). 

As  stated  above,  a  regularization  must  be  used  to  overcome  the  difficulty  associated  with  the 
zero-crossing  in  the  denominator  of  inverse  filter  in  Eq.  (13).  The  zero-crossing  occurs  obvi¬ 
ously  at  a\  =  (k  1  /2) 7T,  k  =  1 , 2,  •  •  • .  The  Tikhonov  regularization  would  replace  1  / cos(ai ) 
in  the  inverse  filter  of  Eq.  (13)  by  the  following 


cos(ai) 
cos(a2)2  +  k 2  ’ 


(14) 


where  K  is  a  small  constant.  However,  this  form  of  regularization  is  not  optimal  since  it  will 
deviate  from  the  original  inverse  filter  when  oci  is  small.  In  other  words,  the  test-object’s  low 
frequency  components  would  not  be  fully  reconstructed  in  the  retrieval.  In  order  to  get  around 
this  drawback,  we  note  that  the  zero-crossing  occurs  only  when  oq  is  close  or  greater  than 
7r/2,  while  the  most  low  frequency  components  have  a  <  1.  Hence  we  set  the  regularization 
parameter  K  =  0  when  a  is  less  than  some  threshold  value  Oq,  and  K  >  0  when  a  is  greater  than 
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Fig.  4.  Comparation  of  the  retrieved  results  using  Eq.  (11)  and  Eq.  (10)  with  Ri  =  5m.  (a) 
is  the  true  Aq;  (b)  is  the  phase-contrast  image  1  acquired  from  Fresnel  diffraction;  (c)  is 
retrieved  using  Eq.  (11);  (d)  is  retrieved  using  Eq.  (10). 


Oq.  In  our  simulation,  we  set  this  threshold  ocq  =  0.8,  so  K  =  0  when  a  <  0.8,  and  jc  =  0.12  when 
a  >  0.8.  We  can  see  from  the  image  in  Fig.  4(c)  that  most  of  the  low  frequency  information 
is  preserved.  On  the  other  hand,  the  quality  of  the  results  can  also  be  numerically  determined 
by  computing  the  error  comparing  to  the  true  Aq  image.  By  taking  the  /^-measure,  which  is 
defined  as  the  integrated  squared  errors,  we  have  the  absolute  errors  for  methods  (10)  and  (11) 
being  0.014  and  0.021  respectively. 

4.  Conclusions 

Phase  retrieval  is  required  for  achieving  artifact-free  x-ray  phase-sensitive  3D  imaging.  Multi¬ 
ple  image-acquisitions  per  projection  are  required  in  general  for  phase  retrieval  [3,  4,  15,  17, 
22,  23,  24].  This  requirement  impedes  the  implement  of  x-ray  phase-sensitive  3D  imaging, 
especially  for  phase-sensitive  3D  imaging  in  clinical  applications.  A  phase-retrieval  approach 
based  on  the  phase-attenuation  duality  with  high  energy  x-rays  can  greatly  facilitate  for  phase 
sensitive  imaging  by  allowing  phase  retrieval  from  only  one  single  projection  image.  However, 
the  previously  derived  phase  retrieval  formula  is  valid  only  for  small  Fresnel  propagator  phases 
corresponding  to  common  clinical  imaging  tasks.  In  this  work  we  generalized  the  previous 
phase  retrieval  formula  to  Eq.  (10),  which  can  be  applied  for  cases  with  large  Fresnel  propaga¬ 
tor  phases  corresponding  to  high  spatial  resolution  imaging.  The  computer  simulation  demon¬ 
strated  superiority  of  this  new  formula  over  the  previous  duality-based  phase-retrieval  formula 
in  reconstructing  the  high  frequency  components  of  imaged  objects.  In  addition,  a  modified 
Tikhonov  regularization  technique  has  been  devised  for  phase  retrieval  in  cases  of  very  high 
resolution  and  large  object-detector  distance  such  that  some  Fresnel  propagator  phases  may  be 
close  or  greater  than  n.  This  new  phase  retrieval  formula  of  Eq.  (10)  lays  the  foundation  for 
implementing  high-resolution  phase-sensitive  3D  imaging  of  soft  tissue  objects. 
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5  paraspinal  SBRT  patients.  Prescription  dose  to  the  PTV  is  24Gy  while 
constraining  the  cord  maximum  dose  to  14Gy.  Plan  evaluation  statistics 
include  PTV  V22.8Gy,  PTV  D95%,  cord  mean  dose,  and  total  beam-on 
time.  Results:  Collimator  rotation  in  the  five  Coll-VMAT  plans  ranges 
from  26°-54°,  with  a  median  of  40°.  Patient-averaged  PTV  V22.8Gy 
(94.6%  with  Coll-VMAT  vs  92.1%  VMAT  and  93.3%  IMRT)  and  D95% 
(22.5Gy  vs  21.4Gy  and  22Gy)  are  higher,  while  cord  mean  dose  is  lower 
(9.8Gy  vs  lO.OGy  and  11.7Gy).  Total  beam-on  time  is  comparable  to 
standard  VMAT  and  substantially  lower  than  IMRT  (5164MU  vs 
4868MU  and  13,283MU).  Conclusion:  The  dosimetric  advantage  of 
Coll-VMAT  is  significant  especially  for  small  targets.  Collimator 
trajectory  optimization  provides  an  additional  degree  of  freedom  in 
VMAT  optimization  that  can  improve  dosimetric  quality  of  paraspinal 
SBRT  compared  to  standard  VMAT  and  IMRT.  Conflict  of  Interest: 
Research  agreement  with  Varian  Corporation 

WE-E-BRB-07 

Quality  Assurance  of  Proton  Range  Compensators  by  CT  Scanning 

Y  Zheng1  *,  E  Klein1  ,  D  Low1  ,  (1)  Washington  University  School  of 
Medicine,  St.  Louis,  MO 

Purpose:  Patient  specific  range  compensators  (RCs)  are  required  in 
passively  scattered  proton  therapy.  Conventionally,  quality  assurance 
(QA)  of  RC  is  performed  by  either  a  visual  inspection  or  spot  check  of 
thicknesses  for  some  drill  points  (e.g.  4).  A  more  comprehensive  QA  of 
RCs  is  needed  for  precise  delivery  of  proton  beams.  The  purpose  of  this 
study  was  to  develop  a  QA  procedure  by  use  of  computed  tomography 
(CT).  Method  and  Materials:  For  the  procedure,  one  or  more  range 
compensators  would  be  scanned  with  a  CT  scanner.  Software  was 
developed  to  automatically  segment  the  range  compensator  on  CT  images 
and  calculate  the  thickness  matrix  for  the  RC.  The  CT-scanned  RC 
thickness  matrix  was  then  compared  to  the  expected  one,  which  was 
reconstructed  from  the  DICOM  ION  Radiotherapy  Plan  exported  from  a 
proton  treatment  planning  system.  The  actual  milling  tool  size  and  shape 
were  taken  into  account  when  the  RC  was  reconstructed.  The  scanned 
RC  thickness  matrix  was  automatically  aligned  to  the  planned  RC  matrix. 
In  addition  to  comparing  the  physical  thickness  of  the  RC,  we  also 
verified  the  water  equivalence  thickness  of  the  range  compensator  based 
on  the  CT  scan.  This  would  allow  us  to  detect  possible  impurities  or 
other  defects  of  the  RC.  As  part  of  the  validation,  we  tested  not  only  an 
unaltered  RC,  but  RC  that  was  damaged  or  altered  post-construction. 
Results:  We  developed  software  to  process  the  scanned  CT  images  of  the 
RC  and  reconstruct  the  planned  RC  from  the  radiation  therapy  treatment 
plan.  We  compared  both  the  physical  thickness  and  water  equivalent 
thickness  between  the  scanned  RC  and  planned  RC  to  quality  assurance 
the  manufacturing  of  the  RC.  Our  process  was  able  to  detect  deviations  of 
less  than  2  mm.  Conclusion:  We  have  developed  a  procedure  to 
efficiently  and  comprehensively  QA  the  RC. 

WE-E-BRB-08 

Motion  and  Delivery  Management  with  Novel  4D-Optimization 

ONohadani1  *,  J  Seco2  ,  T  Bortfeld3  ,  (1)  Massachusetts  General 
Hospital,  Boston,  MA,  (2)  Mass  General  Hospital;  Harvard  Medical, 
Boston  ,  MA,  (3)  Massachusetts  General  Hospital,  Boston,  MA 

Purpose:  Intensity  modulated  radiation  therapy  (IMRT)  has  the  ability  to 
deliver  highly  conformal  dose  distributions  to  tumors  of  complex  shapes. 
However,  motion  uncertainties  can  substantially  degrade  the  quality  of  an 
otherwise  optimized  treatment  plan.  We  present  a  novel  4D-optimization 
technique  that  exploits  time  effects  and  provides  plans  which  are  robust 
against  tumor  motion  and  warrant  deliverability.  Method  and  Materials: 
Our  4D-optimization  technique  takes  all  registered  breathing  phases  into 
account.  All  4D  CT  data  are  conformally  registered.  To  account  for  lung 
tissue  inhomogeneities,  we  use  Monte  Carlo  dose  calculations.  The  4D- 
optimization  method  seeks  to  minimize  the  overall  dose  subject  to  i) 
lower  bounds  for  tumor  voxels  and  ii)  upper  bounds  for  remaining 
voxels,  and  iii)  closeness  of  the  beamlet  intensities  for  neighboring 
phases.  Results:  Based  on  a  clinical  lung  cancer  case,  our  4D- 
optimization  method  delivers  robust  plans  that  significantly  outperform 
the  gating-based  optimal  plans  both  in  terms  of  target  coverage  as  well  as 
OAR  sparing.  In  fact,  the  EUD  to  the  PTV  deviates  only  by  0.1%  from 
the  prescribed  dose  for  our  4D  method,  whereas  the  gated  EUD  deviates 
by  2%.  Moreover,  since  the  4D-robust  plan  is  quasi-static  and  can  be 
delivered  over  all  breathing  phases,  it  warrants  for  a  reduction  in  delivery 
errors  and  is  less  sensitive  to  motion  uncertainties.  Note  that  gating  is 


quite  sensitive  to  the  choice  of  the  gated  phase  and,  thus,  requires 
enlarged  margins.  Furthermore,  due  to  reduced  number  of  segments 
needed,  we  report  an  increase  in  delivery  efficiency.  Conclusion:  Our 
proposed  4D  optimization  method  delivers  plans  that  are  robust  against 
motion  uncertainties.  Moreover,  the  plan  is  also  inherently  robust  against 
sequencing  and  delivery  errors.  Therefore,  it  allows  to  reduce  margins  on 
the  ITV.  The  method  outperforms  gating  both  in  dose  coverage  as  well 
as  delivery  robustness  and  increases  the  delivery  efficiency. 

Joint  Imaging/Therapy  Symposium 
Room  303A 

Image  Guided  Interventions:  Imaging  Guidance 
and  New  Technologies  for  Interventional 
Procedures 

WE-E-303A-01 

Image-Guided  Therapies:  Advances  in  Imaging,  Modeling,  and  New 
Applications 

J  Siewerdsen1  *,  D  Hawkes2  *,  P  Keall3  *,  (1)  University  Health 
Network,  Toronto,  ON,  CA,  (2)  University  College  London,  London, 
GB,  (3)  Stanford  University  Cancer  Center,  Stanford,  CA 

Recent  advances  in  imaging  technology,  image  modeling  and 
registration,  and  advanced  therapeutics  offer  the  potential  for  precise, 
image-guided  therapies  that  will  transform  existing  approaches  to 
planning,  intervention,  adaptation,  and  monitoring  of  therapeutic 
response.  Such  therapies  are  marked  by  a  dramatic  increase  in  the 
application  of  image  information  from  multiple  modalities  and  demand 
accurate  coregistration  of  morphological  and  functional  information 
across  a  broad  range  of  spatial  and  temporal  scales.  The  scope  of 
advanced  therapeutic  modalities  enabled  by  such  advances  is  broad, 
ranging  from  high-precision  image-guided  radiation  therapy  and  surgery 
to  minimally  invasive  target  ablation,  cell-based  therapies,  and  other 
forms  of  novel  therapeutics.  Moreover,  the  information  acquired  during 
the  course  of  intervention  will  drive  patterns  of  therapy  delivery  that  are 
increasingly  adaptive  and  patient-specific.  This  symposium  focuses  on 
the  scientific  principles  and  recent  advances  in  imaging,  guidance,  image 
modeling,  registration,  adaptation,  and  advanced  therapy  delivery 
techniques.  Three  speakers  will  present  on  topics  central  to  such 
advancement:  Dr.  Siewerdsen  will  discuss  principles  of  imaging 
performance,  guidance,  and  the  development  and  integration  of  new 
technologies  for  therapy  guidance;  Dr.  Hawkes  will  address  principles  of 
rigid  and  deformable  image  registration,  fusion,  visualization,  and 
augmented  reality;  finally,  Dr.  Keall  will  discuss  the  latest  advances  in 
imaging,  adaptation,  and  novel  treatment  delivery  techniques  in  image- 
guided  radiation  therapy. 

Learning  Objectives: 

1 .  Gain  exposure  to  new  developments  in  imaging  technologies  for 
image-guided  interventions. 

2.  Understand  challenges  in  rigid  and  deformable  image  registration, 
fusion,  and  visualization. 

3.  Learn  the  latest  advances  in  imaging,  adaptation,  and  therapy 
delivery  techniques  in  IGRT. 

Imaging  Scientific  Session 
Room  304A 

New  Technology,  Techniques  and  Observer 
Performance 

WE-E-304A-01 

Contrast  Mechanism  and  Origin  of  Artifacts  in  Phase-Sensitive  X- 
Ray  Volumetric  Imaging 

X  Wu1  *,  A  Yan1  ,  H  Limj,  (1)  Univ  Alabama  Birmingham,  Birmingham, 
AL,  (2)  Univ  Oklahoma,  Norman,  OK 

Purpose:  Phase-sensitive  x-ray  volumetric  imaging  has  great  potential 
for  tissue  contrast  enhancement  and  radiation  dose  reduction,  compared 
to  conventional  tissue-attenuation  based  x-ray  volumetric  imaging.  In 
order  to  achieve  quantitative  and  artifact-free  phase  sensitive  volumetric 
imaging,  one  should  quantify  the  formed  phase  contrast  and  find  the 
origin  of  artifacts  formed  in  phase-sensitive  volumetric  imaging. 
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Method  and  Materials:  We  incorporated  x-ray  Fresnel  diffraction  into 
x-ray  projection  acquisitions,  and  reconstructed  volumetric  images  from 
these  phase-sensitive  projections  by  the  filtered  backprojection  method. 
A  formula  for  the  reconstructed  tissue’s  apparent  attenuation  coefficients 
was  derived  for  quantitative  analysis.  Computer  simulations  of  phase 
sensitive  3D  tomography  were  performed,  and  reconstructed  3D  phantom 
tomograms  were  analyzed,  and  the  artifacts  associated  with  phase- 
contrast  were  identified.  In  order  to  eliminate  these  artifacts,  we 
performed  phase  retrievals  based  on  the  phase-attenuation  duality  method 
for  each  projection,  and  these  phase-retrieved  projections  were  then  used 
for  new  3D  tomographic  reconstructions.  Results:  The  apparent  linear 
attenuation  coefficients  in  3D  tomograms  were  found  each  being  the  sum 
of  three  parts:  the  tissue  linear  attenuation  coefficient,  the  scaled  3D 
Laplacian  of  tissue  refraction  index,  and  the  artifact  depending  on  global 
distributions  of  tissue  attenuation  coefficients  and  refraction  indices. 
The  simulated  high-resolution  3D  images  of  phantoms  exhibit  severe 
artifacts  as  the  upward  and  downward  density  over-shootings  and  the 
streaks.  The  phase  retrieval  was  shown  very  effective  for  removing  these 
artifacts  and  accurately  reconstructing  tissue  attenuation  coefficients  in 
phase-sensitive  3D  tomography.  Conclusion:  We  clarified  the  contrast 
mechanism  in  phase-sensitive  volumetric  imaging,  found  the  quantitative 
contribution  of  the  phase  contrast  in  reconstructed  3D  images,  and 
revealed  the  origin  of  severe  artifacts.  The  need  of  the  phase-retrievals 
for  quantitative  and  artifact-free  phase  sensitive  volumetric  imaging  was 
clearly  demonstrated. 

Research  supported  in  part  by  DoD  Breast  Cancer  Research  Program 

WE-E-304A-02 

Development  of  Low  Dose  X-Ray  Phase  Contrast  CT  Using  a  Novel 
Projection  Retrieval  Reconstruction  Algorithm 

B  Fahimian1  *,  P  Cloetens2,  Y  Mao1  ,  E  Lee1  ,  J  DeMarco  1  ,  M  McNitt- 
Gray1,  J  Miao1*,  (1)  University  of  California,  Los  Angeles,  CA,  (2) 
European  Synchrotron  Radiation  Facility,  Grenoble,  FR 

Purpose:  In  contrast  to  attenuation  based  x-ray  imaging,  phase  contrast 
x-ray  imaging,  produces  radically  enhanced  images  of  soft  tissues 
through  imaging  of  the  phase  shift  cross-section,  which  is  100-1000 
larger  than  the  corresponding  x-ray  attenuation  cross-section.  One 
obstacle  in  the  development  of  phase  contrast  x-ray  CT  (PCXCT)  for 
medical  imaging  is  the  radiation  dose  imparted  to  the  patient  due  to  large 
number  of  projections  and  flux  required  for  an  adequate  image  quality. 
Here  we  develop  and  quantify  a  methodology  for  low  dose  PCXCT  using 
novel  iterative  Fourier  based  projection  retrieval  reconstruction  algorithm 
that  aims  to  solve  the  reconstruction  problem  from  a  limited  number  of 
projections  through  the  use  physical  and  mathematical  constraints  as  well 
as  the  oversampling  method.  Materials  and  Methods:  An  image  quality 
phantom,  containing  multiple  contrast  visibility  and  resolution  inserts, 
was  imaged  at  the  European  Synchrotron  Radiation  Facility  (ESRF). 
Conventional  reconstructions  were  produced  via  Filtered  Back  Projection 
(FBP)  using  the  full  set  of  200  projections.  Reconstructions  using  the 
projection  retrieval  algorithm  were  made  with  50-60%  fewer  projections. 
The  image  quality  of  the  reconstructions  were  quantitatively  compared 
using  SNR,  CNR,  and  resolution.  Results:  With  60%  fewer  projections, 
the  methodology  matched  the  resolution  of  the  FBP  reconstruction, 
exceeded  the  SNR  of  the  FBP  by  factor  of  1.8,  and  exceeded  the  CNR  of 
the  FBP  by  factor  of  1.7.  No  geometric  deformations  were  observed  due 
to  the  undersampling  of  the  number  of  projections.  Conclusions:  The 
results  indicate  that  the  developed  algorithm  can  reproduce 
reconstructions  of  same  or  superior  quality  using  50-60%  fewer 
projections  and  hence  radiation  dose  than  FBP.  Conflict  of  Interest:  A 
portion  of  this  research  was  sponsored  by  TomoSoft  Technologies,  LLC. 

WE-E-304A-03 

An  Energy-Dispersive  Technique  to  Measure  Tissue  X-Ray  Coherent 
Scattering  Form  Factors 

B  King1  *,  P  Johns1  ,  (1)  Carleton  University,  Ottawa,  ON,  CA 

Purpose:  Data  for  the  x-ray  scattering  properties  of  tissues  are  needed  for 
the  development  of  x-ray  scatter  imaging  systems.  Published  amorphous 
material  form  factors  from  x-ray  diffractometer  measurements  vary 
significantly.  An  energy-dispersive  technique  was  developed  to  provide 
more  accurate  measurements.  Method  and  Materials:  The  form  factor  is 
measured  as  a  function  of  the  momentum  transfer  argument  x=XAsm(6H) 
where  X  is  the  photon  wavelength  and  0  is  the  scattering  angle. 


Crystallographic  diffractometers  measure  scattering  versus  0  at  one 
wavelength  (energy).  Our  method  fixes  0  and  measures  versus  energy, 
using  a  high  purity  germanium  spectrometer  to  measure  the  scattered 
polychromatic  spectrum  from  a  tungsten  anode  x-ray  tube.  The  target  is 
moved  laterally  to  measure  the  transmitted  and  scattered  spectra  without 
adjusting  the  x-ray  tube  and  detector.  By  taking  the  scatter  to  transmitted 
ratio  bin-to-bin,  the  coherent  scattering  form  factor  is  determined  without 
adjustable  parameters.  Repeating  with  four  6  values,  the  range  0.5  nm'1  < 
x  <  12  mu'1  is  obtained.  Results:  Measurements  of  the  coherent  scattering 
form  factor  of  water  and  several  plastics  with  our  prototype  system  give 
good  agreement  for  x  >  7  nm*1  with  the  Independent  Atom  Model.  This, 
together  with  general  agreement  at  smaller  x  with  published 
diffractometer  data,  is  strong  evidence  that  our  measurement  is  accurate. 
The  precision  of  our  measurement  has  been  thoroughly  characterized  and 
is  approximately  10%.  Conclusion:  The  energy-dispersive  technique 
measures  over  a  wider  x  range  than  do  diffractometers  and  does  not 
require  the  use  of  arbitrary  scaling  constants.  Our  measurements  are 
limited  to  5  -  10%  in  x  resolution  at  small  6  but  this  is  acceptable  for 
tissue  measurements.  Resolution  and  precision  will  be  improved  by  using 
rectangular  apertures.  Due  to  the  advantages  of  the  technique  and  to  its 
thorough  characterization,  we  expect  the  energy-dispersive  technique  to 
provide  definitive  form  factor  results. 

WE-E-304A-04 

Local  Determination  of  Tissue  Optical  Properties:  Reconstruction 
Based  On  Fiber-Array  Reflectance  Profiles 

C  Sheng1  *,  C  Bonnerup2  ,  C  Austerlitz1  ,  R  Allison1  ,  C  Sibata1  , 
(1)  Department  of  Radiation  Oncology,  East  Carolina  University, 
Greenville,  NC,  (2)  Department  of  Physics,  East  Carolina  University, 
Greenville,  NC 

Quantitative  evaluation  of  in  vivo  local  tissue  optical  properties  including 
scattering  coefficient  (ps),  absorption  coefficient  (pa)  and  anisotropy  (g) 
is  often  important  in  both  photodiagnosis  and  phototherapy.  In  this  study, 
a  reflectance  based  fast  technique  was  developed  to  determine  the  optical 
properties  of  turbid  media  using  a  linear- array  fiber  bundle  probe.  Five 
200  um  collection  fibers  were  linearly  set  along  from  the  200  um 
illumination  fiber  with  center-to-center  separation  of  350  um.  Spatial 
reflectance  values  were  sequentially  measured  by  spectrometer  connected 
to  a  fiber-switch.  A  model  that  relates  the  reflectance  profiles  to  optical 
properties  of  a  turbid  medium  was  developed  based  on  Monte  Carlo 
simulations  and  phantom  experiments.  Simulation  results  at  wavelength 
of  633  nm  showed  that  ps’  (2~40  cm-1)  and  pa  (0~5  cm-1)  can  be 
determined  by  reflectance  spatial  profiles.  Intralipid  and  Nigrosin  were 
used  to  simulate  different  reduced  scattering  coefficient  (pf)  and 
absorption  coefficient  (pa)  values  within  the  same  range  as  Monte  Carlo 
simulation.  Preliminary  results  show  good  correlation  between  known 
optical  properties  in  tissue  phantom  and  the  measured  optical  properties, 
the  average  error  for  ps’  and  pa  was  7.8%  and  6.6%,  respectively.  With 
same  reduced  scattering  coefficient  (ps’),  changes  in  the  absorption 
coefficient  (pa)  could  be  measured  within  0.1  cm-1.  Accurate  extraction 
of  tissue  optical  properties  from  in  vivo  measurements  could  have 
potential  application  in  noninvasively  superficial  (pre)cancer  detection 
and  phototherapy  planning. 

WE-E-304A-05 

Acoustic  Droplet  Vaporization  for  Enhancement  of  High  Intensity 
Focused  Ultrasound  Thermal  Ablation 

M  Zhang*,  KJ  Haworth,  SD  Swanson,  ML  Fabiilli,  OD  Kripfgans, 
PL  Carson,  JB  Fowlkes,  University  of  Michigan,  Ann  Arbor,  MI 

Purpose:  Acoustic  droplet  vaporization  (ADV)  shows  promise  for  the 
spatial  control  and  acceleration  of  thermal  lesion  production.  In  this 
study,  we  investigated  the  hypothesis  that  ADV  bubbles  can  substantially 
enhance  HIFU  ablation  by  increasing  and  controlling  local  energy 
absorption.  Method  and  Materials:  HIFU/ADV  phantoms  were  made  of 
polyacrylamide,  egg  white,  degassed  water  and  lipid-coated 
perfluoropentane  droplets  (1.8  pm  mean  dia.).  These  egg  white-based 
phantoms  (acoustic  attenuation  of  0.3  dB/cm  at  1.44  MHz)  were  optically 
transparent,  which  allowed  direct  visualization  of  lesion  formation  during 
HIFU  application.  Phantoms  were  exposed  in  degassed  water  at  37°C  to  a 
focused  ultrasound  (1.44  MHz)  with  a  focal  intensity  of  2000  W/cm2. 
Lesion  volumes  were  measured  by  fluid  displacement  or  T2  weighted 
MR  imaging.  Results:  For  2-s  HIFU  exposure  on  droplet-containing 
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